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1 .  INTRODUCTION 

It  is  desired  to  improve  upon  the  ability  to  describe  the  behavior  of 
anisotropic  media  subjected  to  large  pressures,  as  is  the  case  for 
hypervelocity  impact.  It  is  believed  that  expressing  the  anisotropic 
constitutive  relationship  in  a  form  that  makes  use  of  the  deviatoric  stress 
and  strain  tensors  provides  for  a  better  description  of  anisotropic  materials 
whose  compressibility  is  permitted  to  vary  with  volumetric  strain.  The 
deviatoric  stress  technique  is  used  routinely  in  many  impact  codes  for 
describing  isotropic  behavior1-3,  and  is  described  in  many  books  on 
elasticity  and  plasticity^-5.  Anisotropic  schemes  have  also  been 
developed  for  various  impact  codes6-7  which  calculate  a  deviatoric  stress. 
However,  the  deviatoric  stress  is  expressed  in  terms  of  a  total  strain  and  the 
bulk  modulus.  In  a  true  deviatoric  formulation,  deviatoric  stress  is 
expressed  only  in  terms  of  deviatoric  strain,  and  compressibility  affects  only 
the  equation  of  state,  not  the  deviatoric  stress/strain  relation. 

An  anisotropic  formulation  is  proposed  which  satisfies  the  condition  of 
reducing  to  Hooke’s  Law/Prandtl  Reuss  Flow  Rule  when  employing  the  constraint 
of  constant  compressibility  and  isotropy,  but  which  conveniently  allows  for 
anisotropy  and  variable  compressibility.  Additionally,  the  formulation  is 
amenable  for  inclusion  into  existing  impact  codes  which  presently  use  the 
deviatoric  stress  technique  for  isotropic  materials.  A  skeleton  coding  of  the 
scheme  is  provided  in  Appendix  A.  The  scheme  also  provides  an  improved 
technique  for  calculating  hydrostatic  pressure  which  is  less  prone  to  error 
than  existing  techniques.  Finally,  it  is  hoped  that  the  formulation  provides 
an  enhanced  physical  interpretation  on  the  behavior  of  anisotropic  materials 
which  might  otherwise  be  lacking. 


2 .  BACKGROUND 

The  constitutive  relationship  for  any  elastic  material  may  be 
represented  in  contracted  form  as 


where  ct j  and  ej  represent  the  six  independent  stress  and 
strain  components,  and  cij  is  the  modulus  matrix.  The  contracted  form  of 
the  constitutive  relation  is  used  for  the  sake  of  simplicity,  but  the 
tensorial  components  of  the  contracted  form  are  defined  as  follows : 

CTi  =  (ali  a22  ct33  a23  ct13  ct12^ 

€i  :  (£il  e22  e33  e23  e13  CT12> 

In  general,  may  be  a  function  of  a,  e,  e,  etc.  However,  it  is 
somewhat  unwieldy  as  such,  and  is  sometimes  considered  to  be  constructed  of 
constants,  which  produces  the  familiar  Hooke’s  Law.  One  reason  why  the 
deficiency  of  Hooke’s  Law  becomes  apparent  experimentally  under  large 
pressures  is  that  the  bulk  modulus  of  the  material  is  quite  different  from  the 
material’s  stress  free  value. 

For  isotropic  materials,  this  problem  has  been  computationally 
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circumvented  by  the  introduction  of  the  deviatoric  stress  and  strain  tensors. 
These  tensors  differ  from  the  absolute  stress/strain  tensors  in  that  the 
normal  components  of  stress  and  strain  are  decremented  by  the  average  of  the 
normal  stresses  and  strains  respectively.  In  this  way,  the  deviatoric 
quantities  represent  deviation  from  a  hydrostatic  condition,  while  the 
relationship  existing  between  the  average  stress  (negative  of  pressure)  and 
average  strain  (volumetric  dilatation)  is  an  equation  of  state.  Since 
experimental  evidence  reveals  that  the  compressibil ity  of  many  materials 
changes  under  large  pressures,  the  deviatoric  formulation  suggests  that  while 
the  simplicity  of  HooKe’s  Law  (constant  coefficients)  might  possibly  be 
retained  for  computation  of  the  deviatoric  stresses  and  strains,  a  more 
accurate  scalar  equation  of  state  should  simultaneously  be  employed  to  account 
for  non-linear  compressibil ity  effects. 


3.  ELASTIC  DEVIATORIC  ANISOTROPY 


While  the  mathematics  of  the  constant  coefficient  constitutive 
relationship  for  anisotropic  materials  is  well  understood,  the  casting  of 
these  rules  into  a  deviatoric  format  is  not  nearly  as  straightforward  as  it  is 
for  isotropic  materials.  Difficulties  arise  because  of  two  primary 
differences  in  the  behavior  of  anisotropic  materials  with  respect  to  that  of 
isotropic  materials:  (a)  under  hydrostatic  pressure,  strain  is  not  uniform  in 
all  three  directions  of  the  material  coordinates,  and  (b)  except  under 
restrictive  modulus  conditions,  deviatoric  strain  will  produce  volumetric 
dilatation  (i.e.,  two  different  stress  states  with  the  same  pressure  will 
produce  different  dilatations  in  the  material). 


Decomposition  of  the  stress  and  strain  tensors  into  their  hydrostatic  and 
deviatoric  components  yields: 


(2) 


ej 


(3) 


where  are  all  equal  to  the  components  of  hydrostatic  stress 

(a  ;  (at  +  ag  +  03)/3)  for  normal  stress  components  and 

equal  to  zero  for  the  shear  stress  components.  The  term  ej 
represents  the  normal  strains  due  to  hydrostatic  stress,  and  are  formulated  in 
Appendix  C.  One  may  acquire  upon  substitution  into  equation  (1): 


(sA  +  0i  )  =  Cij  (ej  +  e j )  (4) 

where  barred  quantities  represent  conditions  resulting  from  a  hydrostatic 
pressure,  Si  and  ej  are  the  deviatoric  stresses  and  strains 
respectively,  and  is  the  modulus  matrix.  Unlike  the  isotropic 

materials  in  which  a  hydrostatic  pressure  produces  a  uniform  dilatation  in  all 
three  coordinate  directions,  hydrostatic  strain  for  an  anisotropic  material  is 
non-uniform.  Therefore,  if  one  defines  the  deviatoric  components  of  stress 
and  strain  to  be  the  total  stress/strain  components  decremented  by  an  amount 
which  would  result  from  a  hydrostatic  stress  state,  one  can  conclude  (per 
condition  "a"  above)  that  there  is  a  unique  hydrostatic  strain  component 
associated  with  all  three  directions  in  the  material  coordinates  (the 
coordinate  system  which  produces  no  shear  coupling).  Equation  (4)  may  be 
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decoupled  to  give  a  hydrostatic  equation 


and  a  deviatoric  relationship  void  of  hydrostatic  terms: 

si  :  cij  ej  (6) 

For  the  saKe  of  clear  visualization,  the  formulation  will  be  described  for 
transverse  isotropy,  though  extension  to  orthotropy  is  straightforward8. 
Figure  1  depicts  material  elements  from  an  anisotropic  body  whose  material 
(preferred)  coordinate  systems  differ  from  the  laboratory  frame  of  reference. 
The  preferred  coordinate  system  is  the  reference  frame  in  which  the 
constitutive  relation  reduces  to  its  most  simple  form.  Figure  2  shows 
properties  of  the  preferred  transversely  isotropic  material  frame.  Mechanical 
properties  are  invariant  with  respect  to  reference  frame  rotations  that  are 
confined  to  the  plane  of  isotropy.  As  such,  a  certain  symmetry  of  mechanical 
properties  exist  in  transversely  isotropic  materials  which  are  absent  in 
orthotropic  materials.  The  proposed  model  will  be  described  in  the  material 
(preferred)  coordinate  system.  Solutions  of  problems  in  which  the  laboratory 
frame  and  the  material  frame  do  not  coincide  pose  no  problem  if  one  first 
transforms  stress  and  strain  to  the  material  frame  (see  Appendix  B). 


Under  the  influence  of  a  purely  hydrostatic  stress  state  (and  assuming 
the  moduli  to  be  constant),  there  will_be  a  constant  ratio  between  the 
anisotropic  (lon_gitudinal)  strain  and  the  transversely  isotropic 

planar  strain  e2.  Defining  the  ratio  in  terms  of  material 
compliances  sij  (where  sij  =  (Cij)  1): 

ei  Sn  +  2S12 

Ke  =  — -  :  -  (7  ) 

e2  S22  +  S12  +  S23 


it  is  seen  that  this  parameter  (Ke)  reduces  to  a  value  of  unity 
for  isotropy,  where  will  equal  S22,  and  S12  will  equal  S23. 

Using  the  definition  that  deviatoric  stress  is  that  part  of  the  stress 
tensor  which  deviates  from  the  hydrostatic  stress  condition,  one  can  conclude 
that  the  deviatoric  stress  has  no  hydrostatic  component 

Sj  +  s2  +  s3  =  0  (8) 

One  may  substitute  the  deviatoric  constitutive  relation,  equation 
(6) ,  to  acquire 

Ka  e^  +  e2  +  e3  =  0  (9) 

where  Ka  physically  represents  the  ratio  of  longitudinal  and 
transverse  stress  under  conditions  of  uniform  strain  (ej  :  e2  =  £3)1 
and  is  given  by 


cli  +  2C12 

Ka  =  - 

C22  +  C1 2  +  c23 


(10) 
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Figure  1.  The  Preferred  Reference  Frame  of  Material  Elements  May  Not 
With  the  Laboratory  Frame  of  Reference 


Coincide 


4 


4 
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Figure  2.  The  Transversely  Isotropic  Material  Reference  Frame  With  Isotropy  in  the  2-3  Plane. 


As  a  result,  the  sum  of  the  three  normal  deviatoric  strain  increments 
is  not  generally  zero,  but  rather  equals  a  deviatoric  dilatation  (e). 

The  significance  of  this  term  is  that  a  state  of  stress  whose  average  normal 
value  is  zero  can  produce  volumetric  change  on  an  element  with  respect  to  that 
element's  stress  free  volume. 

If  one  wishes  to  convert  a  given  elastic  strain  state  (e^)  into 
the  elastic  deviators  (e^),  elastic  _  deviatoric  dilatation  (e),  and 
the  hydrostatic  strain  components  (ij),  the  following  nine  equations 
given  below  may  be  used  for  a  transversely  isotropic  material  (whose  plane  of 
isotropy  is  the  2-3  plane): 


el  :  el  '  el 

(3a) 

e2  =  e2  "  e2 

(3b) 

e3  =  e3  *  e2 

(3c) 

64  =  £4 

(3d) 

e5  =  e5 

(3e ) 

e6  :  e6 

(3f  ) 

e  =  eA  +  e2  +  e3 

(Dilatation  of  Deviatoric 
Strain ) 

(ID 

®l  :  Ke  e2 

(Non-uniform  hydrostatic 
strain ) 

(7) 

K^ei  +  e2  +  e3  =  0 

(Assures  that  deviatoric 

stress 

has  no  hydrostatic 
component)  (9) 


A  convenient  solution  of  this  set  of  equations  is  given  in  Appendix 
C.  Finally,  the  use  of  the  deviatoric  constitutive  relation,  equation  (6) 
hinged  upon  the  satisfaction  of  equation  (5).  Inverting  equation  (5)  into 
compliance  form  and  summing  the  three  equations  for  normal  strain  yields  upon 
reduction : 

a  :  K  (e^  +  e2  +  e3  -  e )  (12) 


where  K  is  a  true  material  property  which  will  be  called  the  effective 
bulk  modulus  of  the  material  (it  equals  the  reciprocal  of  the  sum  of  the  nine 
normal  compliance  matrix  components),  and  (ej  +  e2  +  63)  is 
the  total  volumetric  dilatation  of  the  material  element.  This  effective 
modulus,  unliKe  the  bulK  modulus,  is  independent  of  deviatoric  stress  in 
anisotropic  materials.  The  bulk  modulus  reduces  to  the  effective  bulk  modulus 
only  when  the  deviatoric  dilatation  e  equals  zero.  This  condition  occurs 
under  either  of  the  following  conditions:  the  material  is  isotropic,  or  the 
loading  is  purely  hydrostatic. 

It  was  mentioned  previously  that  the  empirical  relation  between 
dilatation  and  pressure  is  not  a  linear  one.  One  advantage  of  the  deviatoric 
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formulation  lies  in  the  ability  to  arbitrarily  make  the  hydrostatic  relation 
non-linear  while  retaining  the  linear  simplicity  of  Hooke's  Law  for  the 
deviatoric  portion  of  the  constitutive  relation.  Though  this  ad  hoc  procedure 
does  not  theoretically  follow  as  an  extension  to  Hooke’s  Law,  it  does  permit 
the  code  user  to  more  flexibly  model  the  empirical  behavior  of  the  material. 

There  are  also  codes  employing  the  incremental  strain  approach  which  use 
a  formulation  employing  deviatoric  stress,  though  the  formulation  can  not  be 
termed  deviatoric.  The  form  of  the  relation  used  by  the  HELP  code6  is 


f  C  Ae  -  3K  (Ae  +  Ae  +  Ae  )  ,  i  =  1.  2,  3 
U  j  12  3 

L  cij  AcJ  •  5>  6 


(13) 


where  K  is  identified  as  the  bulk  modulus  which  presumably  can  be  made 
dependent  on  dilatation  (and  therefore  hydrostatic  stress).  In  this  way,  the 
formulation  may  also  provide  the  flexibility  of  a  truly  deviatoric 
formulation.  However,  equation  (13)  is  not  truly  a  deviatoric  relation,  since 
the  deviatoric  stress  increment  is  not  related  to  deviatoric  strain  increment, 
but  rather  is  expressed  in  terms  of  the  total  strain  increment.  The  system  of 
equations  presently  proposed,  equations  (6  and  12),  are  thus  more  attractive 
in  a  theoretical  sense.  Similarly,  it  has  already  been  pointed  out  that  the 
bulk  modulus  (as  opposed  to  the  effective  bulk  modulus  derived  in  equation 
(12))  is  functionally  dependent  on  deviatoric  stress,  and  in  this  sense 
equation  (13)  will  exhibit  flawed  behavior  if  the  deviatoric  variation  in  bulk 
modulus  is  not  modeled.  Finally,  the  flexibility  afforded  in  equation  (13)  by 
allowing  the  bulk  modulus  to  vary  with  hydrostatic  stress  has  the  disturbing 
effect  that  the  resulting  sum  of  the  normal  stress  deviators  is  not  generally 
zero.  If  this  interpretation  of  the  HELP  algorithm  as  described  in  reference 
7  is  correct,  the  use  of  the  term  stress  deviators  to  describe  the  left  hand 
side  of  equation  (13)  does  not  even  seem  justified. 


EPIC7  use  a  form  similar  to  equation  (13)  except  that  K  is  defined 
in  such  a  way  as  to  force  the  sum  of  the  normal  stress  deviators  to  zero. 

This  ad  hoc  procedure  will  coincidentally  mimic  the  behavior  of  equation  (6), 
though  the  formulation  is  in  error  during  the  subsequent  hydrostatic  stress 
calculation  by  not  accounting  for  the  deviatorically  induced  dilatation  (e). 


To  see  additional  advantages  afforded  by  the  proposed  formulation  when 
using  a  code  which  employs  an  incremental  strain  approach,  compare  the 
proposed  algorithm  specifics  with  that  of  the  prior  formulation  used  in 
HELP6.  The  proposed  formulation  takes  strain  increments,  decomposes  them 
into  hydrostatic  and  deviatoric  components.  Equation  (6)  is  used  in  an 
incremental  way  to  update  deviatoric  stress.  If  the  hydrostatic  strain 
increments  are  summed  and  remembered,  equation  (12)  may  be  used  to  evaluate 
the  hydrostatic  stress  value  directly.  If  the  hydrostatic  stress  is  a 
function  of  volumetric  dilatation  only,  then  errors  introduced  into  the 
calculation  of  hydrostatic  stress  are  machine  precision  dependent,  but  not 
algorithm  dependent.  That  is  to  say,  errors  in  the  calculation  of  hydrostatic 
pressure  are  insensitive  to  the  size  of  the  hydrostatic  strain  increment. 


On  the  other  hand,  an  incremental  stress  formulation  like  that  proposed 
for  HELP6  experiences  errors  which  are  dependent  on  hydrostatic  strain 
increment  size  (which  is  proportional  to  the  calculation  timestep  size),  if 
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variable  compressibility  is  employed.  For  example,  use  of  equation  (13)  as 
described  for  materials  with  variable  compressibility  requires  that  some  sort 
of  average  compressibility  be  calculated  for  the  time  increment  in  question. 
As  shown  in  Figure  3,  the  average  bulk  modulus  depends  not  only  on  the  total 
element  dilatation,  but  also  on  the  size  of  the  strain  increment  (since 
dilatation  changes  with  strain  increment).  Therefore,  the  accuracy  of  such  a 
scheme  is  limited  by  the  integration  step  size  regardless  of  machine 
precision.  Presumably,  this  problem  can  be  avoided  if  one  replaces  the 
modulus  dilatation  product  at  the  end  of  equation  (13)  with  a  Aa  term, 
where  the  Aa  term  is  directly  obtainable  Knowing  the  previous  and 
present  cycles’  average  stress. 

However,  many  non-linear  equations  of  state  that  are  routinely  employed 
in  impact  codes  like  HELP&  show  a  dependence  of  hydrostatic  pressure  on 
internal  energy.  Under  such  conditions,  this  dependence  of  pressure  on  energy 
must  effectively  be  reflected  in  equation  (13)  for  consistency  to  be 
maintained.  However,  since  Internal  energy  is  affected  by  the  work  done  by 
the  internal  stresses  (which  include  deviatoric  stresses),  a  coupling  of 
internal  energy,  pressure,  and  deviatoric  stresses  exists.  No  simple  means 
exists  to  solve  this  set  of  equations  simultaneously,  and  a  lengthy  iterative 
process  becomes  necessary.  Since  no  mention  of  such  coupling  and/or  iteration 
was  made  in  reference  6,  it  is  believed  that  none  is  performed.  Thus,  it  can 
be  seen  that  equation  (13)  suffers  many  drawbacks  which  make  its  use  less 
desirable  than  the  proposed  method  given  by  equations  (6  and  12)  in  which  the 
deviatoric  relations  are  free  of  hydrostatic  terms. 

In  summary,  the  steps  proposed  for  deducing  elastic  anisotropic  deviators 
in  equations  (6  and  12)  follow  closely  those  for  isotropic  materials  in  the 
following  ways:  (i)  deviatoric  stress  is  expressible  totally  in  terms  of 
deviatoric  strain,  and  (2)  pressure  is  expressible  totally  in  terms  of 
dilatations. 

The  differences  from  the  isotropic  formulation  may  also  be  noted:  (1)  the 
matrix  relating  deviatoric  stress  to  deviatoric  strain  is  not  diagonal  in  the 
anisotropic  case,  and  (2)  the  total  volumetric  dilatation  must  be  modified  by 
the  deviatorically  induced  dilatation  when  calculating  the  pressure. 


4.  PLASTIC  DEVIATORIC  ANISOTROPY 

The  anisotropic  equivalent  to  the  Prandtl-Reuss  flow  rule  of 
plasticity  can  be  similarly  cast  into  a  deviatoric  form.  Stress  behavior  of 
yielding  material  is  governed  primarily  by  the  nature  of  the  yield  surface, 
which  defines  the  allowable  stress  states  of  the  material  and  subsequent 
plastic  flow  properties  (Appendix  D).  In  general,  only  a  portion  of  a  post¬ 
elastic  strain  increment  (Ae^)  contributes  to  changing  the 
stress.  That  portion  is  designated  the  elastic  strain  increment 
(Ac  j).  The  remaining  portion  of  the  strain  increment  is  designated 
the  plastic  strain  increment  (AejP).  This  decomposition  of  the 
strain  increment  is  governed  by  two  rules:  (1)  an  infinitesimal  plastic  strain 
increment  vector  must  be  normal  to  the  yield  surface  at  the  stress  state  under 
consideration,  and  (2)  a  stress  increment  vector  tending  to  go  outside  of  the 
yield  surface  can  at  most  move  tangentially  to  the  yield  surface  at  the  stress 
state  under  consideration. 
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Figure  3.  Hydrostatic  Pressure  Calculations  Based  on  Dilatation  Increment: 
Have  Their  Accuracy  Limited  by  the  Size  of  the  Dilatation 
Increment. 
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Because  of  the  linearity  of  the  equations  governing  the  conversion  from 
absolute  elastic  strain  (ej)  to  deviatoric  elastic  strain  (e  j, 

e,  e  j),  one  is  assured  that  by  decomposing  the  elastic  strain 
increment  into  any  two  arbitrary  divisions,  the  sum  of  the  two  converted 
strain  divisions  equals  the  conversion  of  the  strain  division  sum.  This  rule 
becomes  handy  for  impact  code  implementation  if  the  two  strain  divisions  are 
taken  as  the  total  strain  increment  and  the  negative  of  the  plastic  strain 
increment  (the  sum  of  which  add  up  to  the  elastic  strain  increment).  In  this 
way,  the  stress  changes  may  be  calculated  on  the  assumption  that  the  total 
stress  increment  is  elastic.  If  it  can  then  be  determined  that  yield  has  been 
violated,  a  fictitious  stress  may  be  calculated  from  the  plastic  strain 
increment,  and  subtracted  from  the  stress  state  which  is  in  violation  of  yield 
to  give  the  true  stress  state. 

To  see  how  this  is  employed  in  actuality,  consider  the  deviatoric 
constitutive  relation,  equation  (6),  in  which  the  deviatoric  stress  increment 
is  calculated  via  the  product  of  the  modulus  and  elastic  deviatoric  strain 
increment.  The  linearity  of  the  deviatoric  conversion  equations  implies,  for 
plastic  deformation,  that: 

Asi  3  Cij  (Ae/  -  Ae/)  (14)  • 

The  deviatoric  total  strain  increment  (Ae  j11)  is  calculated 

with  the  deviatoric  conversion  equations,  based  on  the  total  strain  increment. 
The  plastic  deviatoric  strain  increment  (AejP)  can  be  decomposed 
into  its  total  plastic  (Ae  jP)  and  hydrostatic  plastic 

(AijP)  components  r especti vely. 

The  total  plastic  strain  component  is  necessarily  normal  to  the  yield 
surface,  and  is  given  by: 

3f 

Ae  jP  =  AX  -  (15) 


where  f  is  the  equation  governing  the  yield  surface,  and  Ax  is  a 
proportionality  constant  for  the  yield  surface  normal  (3f/3aj), 
which  has  been  evaluated  at  the  stress  state  in  question.  If  one  assumes  an 
anisotropic  yield  condition  like  Hill’s9  in  which  the  yield  criterion  is 
independent  of  the  hydrostatic  pressure,  then  the  yield  surface  normal  may  be 
evaluated  with  the  use  of  the  deviatoric  stresses  (e.g.  3f/3sj). 

Similarly,  the  hydrostatic  plastic  component  represents  the  three 
components  of  plastic  deviatoric  dilatation,  and  can  be  explicitly  calculated 
knowing  the  elastic  and  plastic  material  constants  and  the  same 
proportionality  constant  AX  required  above. 

As  a  side  note,  the  usage  "plastic  dilatation"  would  seem  to  imply  that 
plastic  incompressibility  does  not  hold.  This  is  however  not  the  case. 

Recall  that  equations  (3,7,9  and  11)  were  proven  valid  only  for  elastic 
deformations.  The  concept  of  plastic  strain  was  introduced  to  represent  the 
difference  between  the  elastic  and  total  strain  components.  This  term 
"plastic  dilatation"  in  fact  represents  a  portion  of  the  total  dilatation  to 
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be  subtracted  off  to  yield  the  proper  value  of  elastic  deviatoric  dilatation. 
The  plastic  incompressibility  relation: 


Ae^P  +  Ae2p  +  Ae3p  =  0 


(16) 


is  still  assumed  to  hold  throughout  all  calculations  derived  here.  Thus, 
expressing  the  plastic  deviatoric  dilatation  term  as 

dejP 

AejP  =  - AA  (17 ) 

dA 


the  deviatoric  constitutive  relation  may  be  expressed,  using  equations 
(14,15,  and  17)  as 


As 


i 


r  m 


V  3s j 


dc  P  ) 
_ j _ 

dA  , 


AX 

j 


(18) 


Notice  that  the  only  term  in  this  relationship  which  differs  from  the 
isotropic  case  is  the  last  term  involving  (dejP/dX).  This 
term  is  zero  for  the  isotropic  case  because  of  the  fact  that  there  is  no 
dilatation  as  a  result  of  deviatoric  stress.  Similarly,  this  term  can  not 
generally  be  zero  for  the  anisotropic  case  because  equation  (18)  is  a 
deviatoric  stress  relationship.  The  term  (dejP/dX)  is 
precisely  the  magnitude  required  to  force  the  deviator  stress  to  remain  in  the 
ir  plane  (i.e.  have  no  hydrostatic  components).  The  derivation  of  this 
term  (dejP/dA)  is  described  in  Appendix  C. 


The  quantity  AX  may  be  evaluated  by  taking  the  scalar  product  of 
equation  (18)  with  (af/asj).  Because  A  s^  is  tangential 

to  the  yield  surface  and  (af/as^)  is  the  yield  surface  normal, 
the  scalar  product  is  zero.  Similarly  the  term  (dejP/dX)  as 

derived  in  Appendix  C  is  of  a  form  identical  to  that  resulting  from  the  purely 
hydrostatic  stress  state_  described  in  equation  (7).  Thus,  it  is  the  case  that 
the  quantity  Cjj(dejP/dX)  is  parallel  with  the  hydrostat 

vector.  If  one  assumes  an  anisotropic  yield  condition  like  Hill’s^  in 
which  the  yield  criterion  is  independent  of  the  hydrostatic  pressure,  the 
scalar  product  of  CijfdejP/dA)  and 
(di/dSy)  is  also  zero.  Thus  the  value  for  AX  may  be 

calculated  as: 


AX 


3f 

asi 


3f 


3f 


as 


(19) 


This  expression  for  AX  is  of  a  form  identical  to  that  obtained  for  the 
isotropic  case,  and  can  be  used  in  equation  (18)  to  calculate  the  elastic 
deviatoric  stress  increment. 
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Because  of  the  curvature  of  the  yield  surface  and  the  fact  that  AX 
was  calculated  for  the  stress  state  existing  at  the  beginning  of  the  time 
cycle,  the  updated  stress  state  resulting  from  equation  (16)  may  in  fact  still 
lie  slightly  outside  the  yield  surface.  What  is  done  at  this  point  in  both 
the  existing  models  and  the  proposed  one  is  to  scale  back  all  the  stress 
components  uniformly  until  the  yield  surface  is  exactly  reached.  Though  this 
technique  introduces  some  error  on  its  own,  it  is  believed  that  the  error  is 
not  too  great  since  the  components  of  the  increment  of  stress  scale  back  are 
nearly  normal  to  the  yield  surface  in  many  cases.  Also,  ways  have  been 
devised  by  Vavrick  and  Johnson7  to  decrease  the  magnitude  of  this  error. 
Their  techniques  employ  subdivision  of  the  time  cycle.  However,  some 
anisotropic  formulations  use  a  deviatoric  stress  formulation  in  which  elastic 
deviatoric  stresses  are  defined  in  the  following  way6 


As 


i 


C  Ae  -  3K  (Ae  +  Ae  +  Ae  )  ,  i  =  i,2,  3 
ij  j  12  3 

,  Cjj  Ae  j  ,  i  =  4,  5,  6 


(13) 


and  additional  error  is  introduced  as  a  result.  This  occurs  because  the 
formulation  in  equation  (13)  does  not  guarantee  that  the  sum  of  the  deviatoric 
stresses  wiil  equal  zero  for  an  anisotropic  material,  and  in  fact  they  will 
generally  not  do  so.  As  a  result,  any  scale  back  of  the  stresses  employed  to 
meet  the  yield  criterion  will  include  a  hydrostatic  component.  Such 
hydrostatic  scale  back  violates  basic  rules  of  yield  surface  normality  in  a 
fundamental  way.  Furthermore,  techniques  proposed  by  Vavrick  and  Johnson 
which  decrease  the  error  resulting  from  stress  scale-back  will  not  decrease 
the  amount  of  hydrostatic  stress  error  introduced  into  the  calculation  as  the 
result  of  using  a  formulation  like  that  of  equation  (13). 


5.  CONCLUSIONS 

An  anisotropic  formulation  has  been  proposed  which  satisfies  the 
condition  of  reducing  to  Hooke’s  Law/Prandtl  Reuss  Flow  Rule  when  employing 
the  constraint  of  constant  compressibility  and  isotropy,  but  which 
conveniently  allows  for  anisotropic  material  properties  and  variable 
compressibility. 

The  deviatoric  stress  technique  which  has  been  used  routinely  in  the 
isotropic  impact  codes  for  describing  isotropic  behavior  has  been  effectively 
combined  with  the  anisotropic  constitutive  relations  to  produce  a  truly 
deviatoric  anisotropic  constitutive  relation.  In  this  deviatoric  formulation, 
deviatoric  stress  is  expressed  only  in  terms  of  deviatoric  strain,  and 
compressibility  does  not  influence  the  deviatoric  relation. 

Existing  formulations  suffer  from  drawbacks  which  have  been  eliminated  in 
the  present  formulation.  Some  of  the  drawbacks  of  previous  formulations  may 
be  enumerated  as  follows:  (1)  working  with  absolute  stress  and  strain  offers 
no  simple  way  to  perform  calculations  involving  variable  compressibility,  (2) 
calculating  hydrostatic  pressure  increments  (instead  of  complete  hydrostatic 
pressure)  can  introduce  error  associated  with  obtaining  and  averaging  the 
tangent  bulk  modulus  over  a  strain  increment  (this  problem  compounded  by  the 
fact  that  Hugoniot  data  is  usually  gathered  in  the  form  pressure  versus 
dilatation,  the  slope  of  which  is  the  tangent  bulk  modulus),  and  (3)  use  of  a 
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"deviatoric"  stress  which  includes  a  hydrostatic  component  will  produce  error 
in  the  pressure  calculation  if  stresses  are  scaled  back  to  satisfy  the  yield 
condition. 

Additionally,  the  formulation  can  be  simply  coded  into  existing  impact 
codes  which  presently  use  the  deviatoric  stress  technique  for  isotropic 
materials.  Finally,  it  is  hoped  that  the  formulation  provides  an  enhanced 
physical  interpretation  on  the  behavior  of  anisotropic  materials. 
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APPENDIX  A 

SKELETON  FORTRAN  CODING  OF  THE  DEVIATORIC  TRANSVERSELY 
ISOTROPIC  ELASTIC  PLASTIC  CONSTITUTIVE  RELATION 
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In  most  explicit  impact  codes,  stress  is  generally  computed  for  a  region 
of  the  mesh  by  providing  a  subroutine  with  the  strain  rates  in  that  region  of 
the  mesh,  the  stresses  in  that  region  of  the  mesh  at  a  previous  time,  and  a 
timestep  over  which  the  strain  rates  act.  Though  each  code’s  constitutive 
relation  routine  use  their  own  unique  notations,  they  all  generally:  (1) 
convert  the  strain  rates  into  deviatoric  strain  rates,  (2)  calculate 
deviatoric  stress  increments  based  on  the  deviator  strain  rates  and  timestep, 
and  increment  the  previous  stress  state  by  this  increment,  (3)  check 
deviatoric  stress  state  for  material  yielding,  (4)  modify  the  deviatoric 
stress  state  to  account  for  plastic  flow  if  necessary,  and  (5)  calculate 
pressure  based  on  the  volumetric  strain,  and  time  increment,  generally  using  a 
non-linear  equation  of  state. 

The  coding  required  to  modify  isotropic  constitutive  subroutines  is 
provided  below,  with  all  variables  defined,  with  hopefully  enough  additional 
comments  to  clarify  where  in  the  old  subroutine  the  new  coding  should  be 
substituted.  The  variable  notations  used  generally  conform  to  those  used  in 
the  EPIC  code2. 
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SUBROUTINE  ASTRES  (REQUIRED  ARGUMENTS) 
c 

c  Anisotropic  stress  increment  formulation: 
c 

c  REVISED  February-June  1985:  Deviator  Anisotropy 
c 

REAL  LAMBDA 
INCLUDE  ’commons. file’ 

COMMON  /ELAST/  SIGK,  EPSK 

COMMON  /YIELD/EPSBAR(1600),BN(3,3),BS(3),CN(3,3),CS(3),MODFLA 
DIMENSION  STR(6),  DFDS(6),  GI(6),  DE(6),  DSIG(6),  SIG(6),  RSG(6), 

&  DEDL(6) 

LAMBDA  =  0. 
c 

c  Generate  required  anisotropic  parameters  if  they  haven’t  been  generated 
c  already, 
c 

IF  (MODFLA  .EQ.  0)  CALL  AGEN 
c 

c  Calculate  anisotropic  deviator  strains  based  on  total  strains 
c 

CALL  DEPS  (I,  ERDOT,  EZDOT,  ETDOT,  EZTDOT,  ERTDOT,  ERZDOT, 

&  DE,  DEPSB) 

c 

c  Compute  rotation  and  change  in  normal  stresses  because  of  rotation 
c 

SPDT  =  SPINRZxDTl 
DSTRN  =  2.xSPDT*SRZ(I) 
c 

c  Obtain  deviator  stresses 
c 

SBAR  =  (SR(I)  +  SZ(I)  +  ST(I))  /  3. 

SRI  =  SR(I)  -  SBAR 

SZ1  =  SZ(I)  -  SBAR 

ST1  =  ST(I)  -  SBAR 

SZT1  =  SZT(I) 

SRT1  =  SRT(I) 

SRZ1  =  SRZ(I) 

C 

c  Strength  variable  SEFF  is  constant  for  my  formulation 
c 

SEFF  =  FU(M) 

C 

c  Transform  stress  to  LTT  frame 
c 

CALL  THETA  (I,  TH) 

CALL  XFORM  (SRI  ,SZ1  ,ST1  ,SZT1  ,SRTi  ,SRZ1, 

&  SIG(1),SIG(2),SIG(3),SIG(4),SIG(5),SIG(6),TH) 

c 

c  Calculate  stress  increment  due  to  element  rotation  (rsg)  in  LTT  frame 
c 

CALL  XFORM  (-DSTRN,  DSTRN,  0.,  0.,  0.,  (SR(I)-SZ(I))#SPDT, 

&  RSG(1),RSG(2),RSG(3),RSG(4),RSG(5),RSG(6),TH) 
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c 

c  Calculate  trial  stress  increment  (dsig)  due  to  strain  changes  (de) 
c  in  LTT  frame 
c 

CALL  CXE  (DE,  DSIG) 
c 

c  Lump  together  strain  induced  stress  (dsig)  and  rotation  induced 
c  stress  (rsg) 
c 

DO  31  K  :  1,  6 

31  DSIG(K)  =  DSIG(K)  +  RSG(K) 
c 

c  Calculate  trial  stress  state 
c 

DO  33  K  :  1,  6 

STR(K)  =  SIG(K)  +  DSIG(K) 

33  CONTINUE 
c 

c  Test  for  yielding 
c 

TERMl  =  0. 

TERM2  =  0. 

DO  35  K  =  1,  3 
DO  34  L  :  1,  3 

34  TERMl  =  TERMl  +  BN(K,L)  *  STR(K)  *  STR(L) 

35  TERM2  =  TERM2  +  BS(K)  *  STR(K  +  3)»*2 
VMISES  =  SQRT(TERMl/2.  +  3.*TERM2) 

IF(VMISES.LE.SEFF)  THEN 

c 

c  Stress  is  elastic.  Transform  stress  bacK  to  RZT  frame... 
c 

SEFF  =  VMISES 
ICHECK(I)  =  0 
DEPSBP  =  0. 

CALL  XFORM  (STR(1),STR(2)1STR(3),STR(4)1STR(5),STR(6), 

&  SR2  ,SZ2  ,ST2  ,SZT2  ,SRT2  ,SRZ2,-TH) 

GO  TO  310 
END  IF 
c 

c  Yield  has  occured:  Determine  ALF,  the  fraction  of  strain 
c  that  is  pre-yield, 
c 

IF  (ICHECK(I)  .EQ.  1)  THEN 
c 

c  Deformation  already  plastic...  elastic  fraction  (alf)  =  0. 
c 

ALF  =  0. 

ELSE 

C 

c  else  must  determine  elastic  fraction  (alf)  (see  VavricK,  Johnson) 
c 

ICHECK(I)  =  1 
TERMl  =  0. 

TERM2  =  0. 
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TERM3  =  0. 

TERM4  =  0. 

TERM5  =  0. 

TERM6  =  0. 

DO  41  K  :  1,  3 
DO  40  L  :  K,  3 
IF  (K  .EQ.  L)  GOTO  40 

TERM1  =  TERM1  -  BN(K,L)  *  (DSIG(K)-DSIG(L))xx2 

TERM3  =  TERM3  -  BN(K,L)  *  (DSIG(K)~DSIG(L))x(SIG(K)-SIG(L)) 

TERM5  =  TERM5  -  BN(K,L)  *  (SIG(K)-SIG(L))xx2 

40  CONTINUE 

TERM2  =  TERM2  +  BS(K)  x  DSIG(K+3)xx2 
TERM4  =  TERM4  +  BS(K)  *  DSIG(K+3)*SIG(K+3) 

TERM6  =  TERM6  +  BS(K)  *  SIG(K+3)xx2 

41  CONTINUE 

AAA  =  TER  Ml/2.  +  3.XTERM2 
BBB  =  TERM3  +  6.XTERM4 
CCC  =  TERM5/2.  +  3.*TERM6  -  SEFFxx2 

ALF  =  (-BBB  +  SQRT(BBBxx2  -  4.xAAAxCCC))  /  (2,  *  AAA) 

END  IF 
c 

c  Calculate  transition  stress  (str)  and  post  elastic  strain  increment  (d 
c 

DO  51  K  =  1,  6 

STR(K)  =  SIG(K)  +  ALFxDSIG(K) 

51  DE(K)  =  (1.  -  ALF)  *  DE(K) 
c 

c  Of  this  post-elastic  strain  increment,  only  that  portion  normal 
c  to  the  yield  surface  is  plastic.  Equation  is 
c 

c  delta  (epsilon  plastic)  =  lambda  x  (df  /d(sigma)) 
c 

c  where  f=constant  functionally  defines  the  yield  surface 
c 

CALL  DFDSIG(STR,  SEFF,  DFDS) 
c 

c  Generate  Cij  (df/d(sigma)j)  vector  (otherwise  known  as  Gi) 
c 

CALL  CXE(DFDS,  GI) 
c 

c  Generate  the  de/dlambda  vector  (based  on  the  transition  stress  str) 
c 

SFACTR  =  STR(l)  /  SEFF 

ETERM  =  -1.5  x  (SIGK-1.)  /  (2.+SIGKxEPSK)  x  BN(1,2)  x  SFACTR 

DEDL(l)  =  ETERM  x  EPSK 

DEDL(2)  =  ETERM 

DEDL(3)  =  ETERM 

DEDL(4)  =  0. 

DEDL(5)  =  0. 

DEDL(6)  =  0. 
c 

c  Calculate  lambda  (happens  to  equal  the  equivalent  plastic  strain) 
c 

TERM1  =  0. 
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TERM2  =  0, 

DO  5<i  K  :  1,  6 

TERM1  =  TERM1  +  GI(K)  x  DE(K) 

54  TERM2  =  TERM2  +  GI(K)  x  DFDS(K) 

LAMBDA  =  TERM1  /  TERM2 
c 

c  Calculate  element  dilitation  resulting  from  plastic  deviator  increment 
c 

DEPSBP  =  -ETERM  *  (2.  +  EPSK)  x  LAMBDA 

c 

c  Since  [lambda  x  (df/do)]  is  the  plastic  strain  vector,  the 
c  elastic  part  of  the  post  elastic  deviator  strain  vector  (LHS)  is  : 
c  [post  elastic  strain  vector  (RHS)]  -  [lambda  x  {(df/do)-(de/dlambda)3] 
c 

DO  56  K  =  1,  6 

56  DE(K)  =  DE(K)  -  LAMBDA  x  (DFDS(K)  -  DEDL(K)) 
c 

c  Multiply  this  elastic  part  of  the  post  elastic  strain  increment  (de) 
c  by  the  modulus  to  find  the  change  in  stress  after  yielding  (dsig) 
c 

CALL  CXE  (DE,  DSIG) 
c 

c  Add  this  actual  stress  change  (dsig)  to  the  transition  stress  (str)  in 
c  order  to  obtain  the  updated  stress  (sig) 
c 

DO  58  K  =  1,  6 

58  SIG(K)  =  STR(K)  +  DSIG(K) 
c 

c  Because  of  the  linear  interpolation  along  the  curved  yield  surface, 
c  a  correction  must  be  made  to  the  stress  to  place  the  stress  back  onto 
c  the  yield  surface 
c 

TERM1  =  0. 

TERM2  =  0. 

DO  60  K  =  1,  3 
DO  59  L  :  1,  3 

59  TERM1  :  TERM1  +  BN(K,L)  x  SIG(K)  x  SIG(L) 

60  TERM2  -  TERM2  +  BS(K)  x  SIG(K+3)x*2 
VMISES  =  SQRT(TERMi/2.  +  3.xTERM2) 

c 

c  Correct  stress  (sig)  to  place  it  back  on  the  yield  surface 
c 

DO  64  K  =  i,  6 

64  SIG(K)  =  SIG(K)  x  SEFF/VMISES 
c 

c  Transform  stress  back  to  RZT  frame 
c 

CALL  XFORM  (SIG(1),SIG(2),SIG(3),SIG(4),SIG(5),SIG(6), 

&  SR2  ,SZ2  ,ST2  ,SZT2  ,SRT2  ,SRZ2,-TH) 

c 

c  EFFECTIVE  PLASTIC  STRAIN 
c 

EBAR(I)  =  EBAR(I)  +  LAMBDA 

c 
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c  Update  dilitation  from  deviator  elastic  and  plastic  calculations 
c 

310  EPSBAR(I)  r  EPSBAR(I)  +  (DEPSB  -  DEPSBP) 
c 

c  modify  dilitation  to  account  for  deviator  stresses  (for  pressure 
c  calculation) 
c 

U  =  U  +  EPSBAR(I) 

c 

c  dW  =  ode 

c  =  ode  +  ode  +  sde  +  se 

c 

c  The  first  two  terms  end  up  in  energy  equation  as  p  dV.  Second  two 
c  terms  appear  below  as  sde. 

c 

SRBAR  =  (SRI  +  SR2) 

SZBAR  =  (SZ1  +  SZ2) 

STBAR  =  (ST1  +  ST2) 

SZTBAR  =  (SZT1  +  SZT2) 

SRTBAR  =  (SRT1  +  SRT2) 

SRZBAR  =  (SRZ1  +  SRZ2) 

FDVMT=DVDOT*DTl/2. 

EDEV  =.5  *  (SRBAR#ERDOT  +  SZBAR#EZDOT  +  STBARxETDOT 
&  +  SZTBARaEZTDOT  +  SRTBARaERTDOT  +  SRZBARaERZDOT) 

&  *  (DVOLI  -  FDVMT  +1.)*DT1 
c 

c  PLASTIC  WORK  FOR  SYSTEM 
c 

IF(ICHECK(I).GT.O)  THEN 

PLAST  =  PLAST  +  (SEFF  *  LAMBDA)*VOL(I) 

END  IF 
c 

c  Calculate  sound  speed  for  eventual  use  in  timestep  calculation 

c 


c 

c  INTERNAL  ENERGY  &  PRESSURE  (use  corrected  dilitation  for  pressure) 
c 


c 

C  NET  STRESSES 
c 


440  SR(I)  : 

SR2  - 

PRES  - 

Q 

SZ(I)  = 

SZ2  - 

PRES  - 

Q 

ST(I)  = 

ST2  - 

PRES  - 

Q 

SRZ(I)  =  SRZ2 
SRT(I)  =  SRT2 
SZT(I)  =  SZT2 
RETURN 
END 
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CXXXXXXXXXXXXXXXXXKXXKXXXKXXXXKXXKXXXXXXXXXKKKXKXXXKXXXXXXXKKXXXXXXXKXKX 

SUBROUTINE  DEPS  (I,  ERDOT,  EZDOT,  ETDOT,  EZTDOT,  ERTDOT,  ERZDOT, 

&  DE,  DEPSBT) 

c 

c  Calculates  the  anisotropic  deviator  strain  increment 
c 

COMMON  /ELAST/  SIGK,  EPSK 
INCLUDE  ’commons. file* 

DIMENSION  DE(6) 
c 

c  Define  6x1  tensorial  total  strain  increment  vector  (in  RZT  frame) 
c 

10  DER  =  ERDOT  xDTl 

DEZ  =  EZDOT  xDTl 

DET  =  ETDOT  xDTl 

DEZT  =  EZTDOT  xDTl  /  2. 

DERT  =  ERTDOT  xDTl  /  2. 

DERZ  =  ERZDOT  xDTl  /  2. 
c 

c  Transform  strain  increment  vector  to  LTT  frame 
c 

CALL  THETA  (I,  TH) 

CALL  XFORM(DER  ,DEZ  ,DET  ,DEZT  ,DERT  ,DERZ, 

&  DE(1),DE(2),DE(3),DE(4),DE(5),DE(6),TH) 

c 

c  Transform  into  deviator  strains,  determine  depsbt  (dilitation  caused  b 
c  total  deviator  strains,  later  to  be  modified  by  plastic  deviatoric 
c  dilitation) 
c 

TERM  =  (SIGK  -  1.)  /  (2.  +  SIGKxEPSK) 

DESUM  =  DE(1)  +  DE(2)  +  DE(3) 

DEPSBT  =  TERM  x  (EPSK  x  DESUM  -  "(2.+EPSK)  x  DE(1)) 

EPST  =  (DESUM  -  DEPSBT)  /  (2.  +  EPSK) 

EPSL  =  EPSK  x  EPST 
DE(1)  =  DE(1)  -  EPSL 
DE(2)  =  DE(2)  -  EPST 
DE(3)  =  DE(3)  -  EPST 
RETURN 
END 

CXKXXXXKKKXKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK 

SUBROUTINE  THETA  (I,  TH) 
c 

c  Calculate  orientation  of  element  by  any  appropriate  means 
c 


RETURN 

END 

CKKKKKKKKXKXKXKKXXKKKXKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKXXKKKKKKKKKKKKKKKKK 

SUBROUTINE  CXE  (EE,  SS) 
c 

c  Multiplies  on  axis  modulus  by  vector  EE  to  obtain  vector  SS 
c 
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COMMON  /YIELD/  EPSBAR(160O),BN(3,3),BS(3),CN(3,3),CS(3),MODFLA 
DIMENSION  EE(6),  SS(6) 

DO  20  I  =  1,  3 
SS(I)  =  0. 

DO  10  J  =  i,  3 

10  SS(I)  =  SS(I)  +  CN(I,J)  *  EE(  J) 

20  SS(I+3)  =  CS(I)  x  EE(I+3) 

RETURN 

END 

CXKXKXKXKXXXXXXKKXXKKXKXKKKXXXKXXKXXKKXXKKKKKXKXKXKKXXXXKKXKXKXKXXKKKXKX 

SUBROUTINE  AGEN 
c 

c  Calculate  the  on-axis  modulus  matrix  and  yield  parameters  once  only 
c 

COMMON  /ELAST/  SIGK,  EPSK 
COMMON  /ORIENT/  ANGLE,  TPARAM(1600) 

COMMON  /YIELD/  EPSBAR(1600),BN(3,3),BS(3),CN(3,3),CS(3),MODFLA 

COMMON/LUS/LUI,LUP,LUT,LUPR,LUST,LUFAST 
DATA  LUA  /13/ 
c 

OPEN  (LUA,  FILE=’amatl.dat’,  STATUS-’old’) 

REWIND  (LUA) 

MODFLA  :  1 
WRITE  (LUP,  505) 

505  FORMAT  (//’  Calculating  Anisotropic  Modulus’//) 
c 

c  Engineering  Constants:  (for  trans. -isotropic  material) 
c 

c  Longitudinal  Young’s  Modulus 
READ  (LUA,  x)  EL 
c  Transverse  Young’s  Modulus 
READ  (LUA,  x)  ET 

c  Shear  Modulus  in  Longitudinal-Transverse  Plane 
READ  (LUA,  x)  GLT 

c  Shear  modulus  in  transverse  (isotropic)  plane 
READ  (LUA,  x)  GTT 
c  Bulk  Modulus: 

READ  (LUA,  x)  FK 

c  Left  to  calculate:  isotropic,  LT,  and  TL  Poisson  Ratios 
VTT  =  ET/(2.xGTT)  -  1. 

VLT  =  .25  +  (l.-VTT)xEL/(2.xET)  -  EL/(4.xFK) 

VTL  =  VLT  x  (ET/EL) 
c 

c  Modulus  Matrix  (transversely  isotropic): 
c 

DEL  =  (1  -  2xVLTxVTL  -  VTTxx2  -  2xVLTxVTLxVTT)  /  (EL  x  ETxx2) 

CL  =  (1  -  VTTxx2  )  /  (ETxx2  x  DEL) 

CT  =  (1  -  VTLxVLT)  /  (EL  x  ET  x  DEL) 

CLT  =  (VLT  +  VTTxVLT)  /  (EL  x  ET  x  DEL) 

CTT  =  (VTT  +  VLTxVTL)  /  (EL  x  ET  x  DEL) 

CG  =  GLT 
C 

SL  =  1.  /  EL 
ST  =  1.  /  ET 
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SLT  =  -VLT  /  EL 
STT  =  -VTT  /  ET 
SG  =  1.  /  (2.  n  GLT) 
SGI  =  ST  -  STT 


c 

c  Calculate  Keps  and  Ksig  (variables  EPSK  and  SIGK  respectively) 
c 

EPSK  =  (SL  +  2.#SLT)  /  (ST  +  SLT  +  STT) 

SIGK  =  (CL  +  2.xCLT)  /  (CT  +  CLT  +  CTT) 

c 

CN(1,1)  =  CL 
CN(1,2)  =  CLT 
CN(1,3)  =  CLT 
CN(2,2)  =  CT 
CN(2,3)  =  CTT 
CN(3,3)  =  CT 
CS(i)  =  (CT  -  CTT) 

CS(2)  :  2,  »  CG 

CS(3)  =  2.  «  CG 

CN(2,1)  =  CN(1,2) 

CN(3,2)  =  CN(2,3) 

CN(3,1)  =  CN(i,3) 

WRITE  (LUP,  x)  ’Compliance  Matrix:’ 

WRITE  (LUP,  10)  SL,  SLT,  SLT 

WRITE  (LUP,  10)  SLT,  ST,  STT 

WRITE  (LUP,  10)  SLT,  STT,  ST 
WRITE  (LUP,  10)  SGI,  SG,  SG 

WRITE  (LUP,  11) 

WRITE  (LUP,  X)  ’Modulus  Matrix:  ’ 

WRITE  (LUP,  10)  CN(1,1),  CN(1,2),  CN(1,3) 

WRITE  (LUP,  10)  CN(2,1),  CN(2,2),  CN(2,3) 

WRITE  (LUP,  10)  CN(3,1),  CN(3,2),  CN(3,3) 

WRITE  (LUP,  10)  CS(i),  CS(2),  CS(3) 

10  FORMAT  (3(E15.7,4X)) 

C 

c  Read  Orientation  of  anisotropy 
c 

READ  (LUA,  x)  ANGLE 
c 

c  Calculate  Yield  parameters  (bn(i,j)  ,  bs(i))  (transversely  isotropic) 
c 

c  Longitudinal  Strength 
READ  (LUA,  x)  SIGL 
c  Transverse  Strength 

READ  (LUA,  x)  SIGT 
c  LT  Shear  Strength 

READ  (LUA,  x)  SIGLT 
c 

SEFF  =  SIGT 

BN(1,1)  =  2.  x  SEFFxx2  /  SIGLxx2 
BN(2,2)  :  2.  x  SEFFxx2  /  SIGTxx2 
BN(3,3)  =  BN(2,2) 

c 

BN(i,2)  =  -(+BN(1,1)  +  BN(2,2)  -  BN(3,3))  /  2. 
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BN(1,3)  =  -(+BN(1,1)  -  BN(2,2)  +  BN(3,3))  /  2. 

BN(2,3)  =  -(-BN(l.l)  +  BN(2,2)  +  BN(3,3))  /  2. 

BN(2,i)  =  BN(1,2) 

BN(3,1)  =  BN(1,3) 

BN(3,2)  =  BN(2,3) 

c 

FACTOR  =  SIGL  /  SIGT 

TAU2  =  SIGT##2  /  (4.  -  (l./FACTORxx2)) 

BS(i)  =  SEFFxx2  /  (3.  *  TAU2) 

BS(2)  =  SEFFxx2  /  (3.  x  SIGLTxx2) 

BS(3)  =  BS(2) 

WRITE  (LUP,  ii) 

WRITE  (LUP,  x)  ’Strength  Matrix:’ 

WRITE  (LUP,  10)  BN(i,l),  BN(1,2),  BN(1,3) 

WRITE  (LUP,  10)  BN(2,1),  BN(2,2),  BN(2,3) 

WRITE  (LUP,  10)  BN(3,1),  BN(3,2),  BN(3,3) 

WRITE  (LUP,  10)  BS(1),  BS(2),  BS(3) 

WRITE  (LUP,  11) 

11  FORMAT  (/) 

RETURN 

END 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

SUBROUTINE  XFORM  (U1,U2,U3,U4,U5,U6,P1,P2,P3,P4,P5,P6,TH) 
c 

c  Transforms  stresses  and  strains: 

c  ui  :  stress  or  strain  prior  to  transformation  (unprimed  frame) 

c  pi  :  stress  or  strain  after  transformation  (primed  frame) 

c  th  :  CCW  angle  of  transformation  (in  RZ  frame) 
c 

COMMON/LUS/LUI, LUP, LUT.LUPR, LUST, LUF  AST 
c 

REAL  M,  N,  M2,  N2,  MN 
C 

M  =  DCOS(TH) 

N  =  DSIN(TH) 

M2  :  Mxx2 

N2  =  Nxx2 

MN  =  MxN 
C 

c  All  transformations  are  tensorial,  so  stress  and  strain  are  the  same 
c 

PI  =  M2xUl  +  N2xU2 

P2  =  N2xUl  +  M2xU2 

P3  = 

P4  = 

P5  = 

P6  =  -(MN)xUl  +  (MN)xU2 

C 

RETURN 
END 

CXKXXXKXXXKXXXXXXXXXKXXXXXXXXXXKXXXXXXXXXXXXXXXXKXKXXXXXXXXKXXXXXXXXKXXX 

SUBROUTINE  DFDSIG  (S,  SEFF,  DFDS) 
c 

c  Calculates  df/d(sigma)  for  element  in  question 


+  (2.xMN)xU6 
-  (2.xMN)xU6 
U3 

MxU4  -  NxU5 
NxU4  +  MxU5 

+  (M2-N2)xU6 
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COMMON  /YIELD/  EPSBAR(1600),BN(3,3),BS(3),CN(3,3).CS(3),M0DFLA 
DIMENSION  S(6),  DFDS(6) 

TWOS  =  2.*SEFF 


DFDS(l)  = 

(-BN(1,2)»(S(1)-S(2))  -  BN(1,3)*(S(1)-S(3))) 

/ 

TWOS 

DFDS(2)  = 

(  BN(1,2)«(S(1)-S(2))  -  BN(2,3)*(S(2)-S(3))) 

/ 

TWOS 

DFDS(3)  = 

(  BN(1,3)«(S(1)-S(3))  +  BN(2,3)#(S(2)-S(3))) 

/ 

TWOS 

DFDS(4)  = 

3.#BS(1)*S(4)  /  TWOS 

DFDS(5)  = 

3.*BS(2)*S(5)  /  TWOS 

DFDS(6)  = 

RETURN 

END 

3.#BS(3)#S(6)  /  TWOS 

29 


APPENDIX  B 

THE  EFFECT  OF  MATERIAL  FRAME  ON  ANISOTROPIC  COMPUTATIONS 
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Material  frame  is  not  a  consideration  in  isotropic  codes,  because  the 
constitutive  relation  is  identical  in  all  reference  frames.  As  such,  doing 
the  calculations  in  the  laboratory  frame  of  reference  is  the  logical  choice. 
However,  when  anisotropy  is  involved,  the  material  properties  are  different  in 
different  reference  frames.  For  regular  types  of  anisotropy  (e.g.  transverse 
isotropy,  orthotropy,  etc.),  there  are  preferred  directions  in  which  the 
materials  constitutive  relations  reduce  to  their  most  simple  forms.  In 
general,  this  material  frame  does  not  coincide  with  the  laboratory  frame  of 
reference.  Unfortunately,  it  is  usually  the  laboratory  frame  in  which  system 
properties  (stress,  strain,  etc.)  are  described.  Two  approaches  may  thus  be 
taken  to  implement  anisotropy  into  the  codes:  i)  transform  laboratory  stress 
and  strain  into  the  material  frame,  perform  constitutive  computations  in  the 
material  frame,  and  transform  the  resulting  stresses  and  strains  back  into  the 
laboratory  frame,  or  2)  transform  the  simple  material  frame  constitutive 
relations  into  the  laboratory  frame  of  reference,  and  perform  calculations 
with  these  new  laboratory  frame  constitutive  relations. 

The  following  is  a  comparison  of  the  pertinent  relations  as  they  would 
appear  in  both  the  material  and  laboratory  frame  coordinate  systems.  In  Table 
B-i,  primed  values  of  stress  and  strain  denote  values  in  the  laboratory  frame, 
while  unprimed  values  denote  the  material  frame  values.  The  relationship 
between  material  and  laboratory  frame  stress  and  strain  is 


°1  =  TU  ’ 


e  =  T  e  .  > 

1 


(B-i) 

(B-2) 


where  TAj  is  the  appropriate  transformation  matrix  between  laboratory 
and  material  coordinate  systems.  Note  that  because  the  contracted  stress  and 
strain  notations  are  being  used,  the  transformation  matrix  Ty  is  not 
symmetric. 


Table  B-l  shows  the  nature  of  the  calculations  when  done  in  both  the 
laboratory  and  material  reference  frames.  In  the  Table  B-i,  the  substitution: 


3f 


has  been  made  for  simplicity  of  transformation,  where  f  is  the  function 
defining  the  yield  surface  and  3f/3ai  is  the  vector  normal  to 
the  yield  surface  in  the  material  reference  frame.  Table  B-2  contains  the 
general  form  of  the  terms  contained  in  Table  B-l.  For  calculations  done  in 
the  material  frame,  there  is  a  constant  "overhead"  penalty  of  making  the 
initial  stress  and  strain  transformations,  which  does  not  exist  in  the 
laboratory  frame  scenerio.  However,  it  can  be  seen  from  Table  B-i  that  in  the 
laboratory  frame  there  is  penalty  of  transformation  for  every  calculation 
done. 


Thus,  if  anything  but  the  most  trivial  of  calculations  are  required,  then 
it  computationally  pays  to  first  transform  stress  and  strain  to  the  material 
frame,  perform  the  calculations  there,  and  transform  back  at  the  conclusion  of 
the  computations. 
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Table  B-l.  Comparison  of  Governing  Equations  in  the 
Material  and  Laboratory  Coordinate  Frames 


Transformation  to 
Desired  Frame 


Elastic  Constitutive 
Equation 


Yield  Equation 


Plastic  Strain 


Plastic  Strain  AX 

Parameter 


Material  Frame 
(a)  =  [T]{a»l 
(el  =  [T] { e ’ } 

{a)  =  [C]  (e) 

f2  =  (a) 

(AeP)  =  AX 03 {al 

fa)T[<p]T[C]  {Ael 

-  AX 

(cr)T[<i>]TtC]  [<p]  {al 


Laboratory  Frame 
N/A 
N/A 

{a'l  =  [T]-1[C]  [T]{e>) 

f2  =  {a’  )t[T]t[9]  [T]  la’  1 

{Ae'P}  =  AX[T]_1[<p]  [T]  {a>) 

{a’  )T[T]T[(p]T[C]  [T]  {Ae 1 1 
{a’  1 T [T] T [«f»3 T [C]  [<p]  [THa’l 


Transformation  to  {a’l  =  [T]_1{a)  N/A 

Original  Frame 


where: 

[  ]  denotes  a  6x6  matrix, 

{  )  denotes  a  6x1  vector, 

the  superscript  -1  denotes  a  matrix  inverse, 

the  superscript  T  denotes  a  matrix  transpose,  and 

the  vectors  and  matrices  used  in  this  table  are  defined  in  Table  B-2. 
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Table  B-2.  General  Forms  of  Pertinent  Orthotropic  Terms 


Constitutive  Matrix  [C] 


f  c 

11 

C 

12 

C 

13 

0 

0 

0 

c 

12 

C 

22 

C 

23 

0 

0 

0 

C 

13 

C 

23 

C 

33 

0 

0 

0 

0 

0 

0 

C 

44 

0 

0 

0 

0 

0 

0 

C 

55 

0 

,  o 

0 

0 

0 

0 

C66 

Yield  Normal  Matrix  |>] 


i/2 


(  B  B  B  0 


11 

12 

13 

B 

12 

B 

22 

B 

23 

0 

B 

13 

B 

23 

B 

33 

0 

0 

0 

0 

3B 

44 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

3B66  - 


where: 

Bi2  =  (~Bli  "  b22  +  b33  )  /  2 

B13  =  (_Bii  +  b22  -  b33  )  /  2 

b23  -  (  B11  "  b22  ~  b33 )  /  2 


Stress  Vector 


fa]T 


■> 

o  a  o  c  c  a 
11  22  33  23  13  12 


Strain  Vector 


■> 

e  e  e  e  e 
22  33  23  13  12 
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APPENDIX  C 

DERIVATION  OF  GOVERNING  RELATIONS 
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I.  Elastic  Strain  Decomposition: 


For  transversely  isotropic  material,  with  isotropy  in  the  2-3  plane,  the 
decomposition  of  a  given  elastic  strain  state  (e^)  into  the  elastic 
deviatoric  strains  (e^,  elastic  deviatoric  dilatation  (e),  and  the 
hydrostatic  strain  components  (c  j)  has  been  shown  to  require  the 
solution  of  the  following  nine  equations : 


(3a) 
(3b) 
(3c) 
(3d) 
(3e ) 
(3f  ) 

(Dilatation  of  Deviatoric 
Strain  (11) 


ei  :  Ke  e2 


(Non-uniform  Hydrostatic 
Strain)  (7) 


Ka  ei  +  e2  +  e3  =  ° 


(Assures  that  deviatoric 
stress  has  no  hydrostatic 
component)  (9) 


Standard  equation  reduction  techniques  may  be  employed  to  obtain  the 
following  solution  sequence: 


(K, 


e  = 


(2  +  K 


el  + 


-2  ' 


el  :  Ke  e2 


Equations  (3) 


-  1) 


aKe> 


f 


v. 


e2  +  e3  ~ 
(2  +  Ke) 

are  now 


1 


(ej  +  e2  +  e3)  -  (2  +  K£)  €jJ 

Kq  et  ♦  eg  ♦  e3 
(2  +  Kq  Ke) 

directly  solvable  for  ej. 


(C-l) 


(C-2) 

(C-3) 
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II.  Derivation  of  (dejP/dX): 

In  determining  the  elastic  strain  components  to  be  used  for  the 
calculation  of  stress  after  yielding,  it  was  found  to  be  convenient  to 
decompose  these  elastic  components  into  the  total  strain  increment,  and  the 
negative  of  the  plastic  strain  increment.  The  plastic  flow  relations  required 
the  knowledge  of  the  term  AejP,  which  at  the  time  was  left  only 
as  (dejP/dX)AX,  the  quantity  AX  being  determined 

through  "  other  means.  The  term  (dejp/dX)is  acquired  by 
employing  the  deviatoric  conversion  equations  (C-l.C-2,  C-3,  and  3)  on  the 
plastic  portion  of  the  strain  increment.  Again,  this  is  permitted  because  of 
the  linearity  of  the  conversion  equations,  the  negative  of  the  plastic  strain 
increment  being  nothing  more  than  a  decomposed  component  of  the  elastic  strain 
increment. 

Employing  equation  (C-i)  and  making  use  of  the  plastic 
incompressibility  relation  (16),  the  dilative  quantity  AeP  is 
determined  to  be : 

-(Ka  -  1 )(2  +  Ke) 

AeP  =  - Ae^  (C-4) 

(2  +  Ka  Ke> 


Employing  the  first  order  approximation  to  the  plastic  flow  rule,  one 
acquires  AejP  =  AXOf/ao^,  Using  the  relations 
of  Appendix  D  under  the  constraints  of  a  transversely  isotropic  material,  one 
can  show 


3f  -3  Si 


2 


(C-5) 


Thus,  for  the  transversely  isotropic  material  in  question,  the  dilative 
quantity  Aep  may  be  cast  completely  in  terms  of  available 
quantities  (excepting  AX)  as: 


3(Ka  -  i )  (2  +  Kg) 

AeP  =  -  Big  Si  AX  (C-6) 

2(2  +  Ka  Ke ) 

Equations  (C-2.C-3  and  C-6)  may  then  be  employed  to  ascertain  the  quantity 
AejP  as 


-3(K  -  1)  B  s 

Ae  P  =  - ^ - 12—i 

J  2  (2  +  K  K  ) 

a  e 


f  K  > 
e 


L  i  J 


AX 


(C-7) 


As  a  result  of  differentiating  equation  (C-7),  the  quantity  (dejP/dX) 
is  readily  available  for  use  in  equations  (17)  and  (16). 
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APPENDIX  D 

ANISOTROPIC  YIELD  AND  FLOW  RULE  RELATIONS 
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The  theory  and  computer  code  implementation  of  yield  and  plasticity  rules 
for  anisotropic  materials  has  been  detailed  by  others&'7’’.  The  theory 
extends,  the  approach  of  the  Von  Mises  yield  criterion,  which  is  used 
extensively  for  isotropic  materials.  A  simple  review  of  the  pertinent  points 
will  be  done  just  for  clarity.  Hill’s  original  statement  of  the  anisotropic 
yield  criterion  was  given  as : 

2F  =  1  =  E(a2~a3)2  +  GfOj-Cj)2  +  Hfaj-Og)2  + 

2(La42  +  Ma52  +  Na&2 )  (D-l ) 

By  making  the  appropriate  substitutions,  this  criterion  was  restated  by 
Vavrick  and  Johnson  as: 


f  =  i 


•5(Bllal2  + 


B22a22  + 


B33a32)  + 


B12ala2  + 


B13ala3  + 


B23a2a3 


3 


B55a52  + 


B66a62  > 


1/2 


(D-2) 


The  yield  function  f,  when  fixed  at  a  value  of  unity,  implies  a  perfectly 
plastic  material.  Uniform  work  hardening  may  be  realized  by  letting  the  yield 
function  f  take  on  values  greater  than  unity.  The  form  of  equation  (D-2) 
makes  it  easy  to  define  the  material  constants  of  the  B  matrix.  If  for 
example,  Y ^  represents  the  tensile  strength  of  the  material  in  material 
direction  1,  then  considering  the  simple  case  of  uniaxial  tension  in  the  1 
direction,  substitution  into  (D-2)  reveals  directly  that: 

Bn  :2f2/  YA2  (D-3) 


Similarly,  the  other  constants  are  generated  easily  from  simple 
tension  and  shear  data,  or  from  linear  combinations  of  the  other  constants. 
For  transversely  isotropic  materials  such  as  the  ones  being  considered  in  this 
report ,  the  yield  matrix  Bi j  takes  the  form : 


rBil 

B12 

B12 

0 

0 

0  1 

b12 

B22 

B23 

0 

0 

0 

B12 

B23 

B22 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

B55 

0 

,  o 

0 

0 

0 

0 

B55  j 

where: 

Bn  /  V 

B22  =  2  f2  /  Y22 


(D-4 ) 
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Bi2  =  _Bil  /  2 


®23  =  (®li  ~  2B22  )  /  2 
B44  =  f2  /  3  Y42 
B55  --  f2  /  3  Y52 

Because  of  the  fact  that  the  2-3  plane  is  isotropic,  the  yield  criterion 
must  be  independent  of  rotations  in  that  plane.  By  evaluating  the  yield 
equation  (D-2)  under  conditions  of  pure  shear  in  the  2-3  plane,  and  by  then 
reevaluating  yield  in  a  coordinate  frame  rotated  by  45  degrees  in  the 
isotropic  plane,  it  can  be  shown  that  Y4  is  constrained  for  transversely 
isotropic  materials  to  be : 


Y42 


4  - 


2 


t  Yi  J 


(D-5) 


For  calculations  involving  the  yield  surface  normal,  partial  derivatives 
are  taken  on  equation  (D-2)  with  respect  to  each  of  the  nine  tensorial  stress 
components,  and  evaluated  at  the  stress  state  in  question.  Because  the  yield 
equation  (D-2)  is  expressed  in  terms  of  a  convenient  (albeit  non-tensorial) 
six  dimensional  contracted  "stress  vector"  space,  it  must  be  realized  that 
terms  like  04  in  reality  represent  the  sum  of  two  equal  valued  shear 
stresses  (e.g.  =  .5(a23  +  ct32)).  As  such  the 

partial  derivatives  with  respect  to  the  shear  stresses  is  one  half  that  if 
calculated  strictly  on  the  basis  of  equation  (D-2). 
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SYMBOLS 


(  )’  a  primed  variable  is  a  quantity  whose  value  is  taken  in  an  arbitrary 

laboratory  reference  frame.  Unprimed  quantities  are  those  taken  in  the 
"material  coordinate  frame"  of  a  transversely  isotropic  material. 

(  )*  a  superscript  t  denotes  that  a  variable  represents  a  total  quantity, 
which  is  composed  of  an  elastic  part  and  a  plastic  part. 

(  )P  a  superscript  p  denotes  that  a  variable  represents  a  plastic  quantity. 

A (  )  a  delta  before  a  quantity  signifies  that  the  quantity  is  an  increment. 
As  with  all  models  employing  the  Cauchy  strain  tensor,  incremental 
constitutive  relations  must  be  employed  with  corrections  for  rotation 
in  order  to  make  the  proposed  model  acceptable  for  computation  of 
systems  involving  large  strains. 

C’ij  modulus  matrix  (6x6)  which  relates  stress  components  ai  ’  to  strain 
components  ej ’ .  The  "material  coordinate  frame"  of  a  transversely 
isotropic  material  will  be  defined  as  the  reference  frame  whose  the 
modulus  matrix  (designated  without  the  use  of  primes)  is: 

fCCCOOO'l 


11 

12 

12 

c 

12 

C 

22 

C 

23 

0 

0 

0 

c  = 

C 

12 

C 

23 

C 

22 

0 

0 

0 

ij 

0 

0 

0 

C 

44 

0 

0 

0 

0 

0 

0 

C 

55 

0 

,  o 

0 

0 

0 

0 

C 

ctx  elastic  stress  components  in  contracted  notation;  indices  1  to  3  are 
normal  components,  whereas  4  to  6  are  the  shear  components  23,  13  and 
12  respectively. 

e j  elastic  strain  components  in  contracted  notation;  indices  1  to  3  are 

normal  components,  whereas  4  to  6  are  the  shear  components  23,  13  and 
12  respectively . 

Oi  average  stress,  by  definition  equal  to  the  negative  of  the  hydrostatic 
pressure. 

sx  deviatoric  elastic  stress  components  (6  independent).  In  this  report, 
the  term  "deviatoric"  will  imply  a  deviation  from  the  stress  state 
resulting  from  a  condition  of  hydrostatic  pressure. 

ej  deviatoric  elastic  strain  components  (6  independent).  In  this  report, 
the  term  "deviatoric"  will  imply  a  deviation  from  the  strain  state 
resulting  from  a  condition  of  hydrostatic  pressure.  For  anisotropic 
materials,  strain  is  not  uniform  under  conditions  of  hydrostatic 
pressure  (i.e.  the  three  principal  components  of  strain  are  not 
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identical).  As  a  result,  the  normal  deviatoric  strain  components  are 
NOT  simply  the  difference  between  the  total  strain  component  and  the 
average  of  the  normal  strain  components. 

e  deviatoric  dilatation  (e1+e2+e3).  Though  dilatation  is  only  a 

function  of  pressure  for  isotropic  materials,  dilatation  may  vary  in  an 
anisotropic  material  just  by  varying  the  deviatoric  stress  (without 
changing  the  pressure).  Thus,  this  dilatation  associated  with  the 
deviatoric  stress  is  is  referred  to  as  deviatoric  dilatation. 

ij  strain  state  resulting  from  hydrostatic  pressure.  For  an  isotropic 

material,  the  three  normal  "hydrostatic"  strains  would  be  equal.  This 
is  not  the  case  for  anisotropic  material. 

Ka  a  parameter  which  represents  the  ratio  of  longitudinal  to  transverse 
strain  (in  the  material  reference  frame)  under  conditions  of 
hydrostatic  pressure  (c^  =  og  =  <73). 

Ke  a  parameter  which  represents  the  ratio  of  longitudinal  to  transverse 
stress  (in  the  material  reference  frame)  under  conditions  of  uniform 
strain  (ej  =  e2  =  e3). 

9f/3aj  the  vector  normal  to  the  yield  surface,  which  is  given  by  the  function 
f. 

3f/9sj  is  equivalent  to  3f/3cj  for  a  yield  criterion  like  the  Von  Mises 
or  Hill,  where  yielding  is  not  a  function  of  hydrostatic  pressure. 

AX  a  proportionality  constant  between  the  yield  surface  normal  vector,  and 
the  total  plastic  strain  increment  vector,  which  are  parallel. 

axial  flow  stress  along  the  longitudinal  material  direction  (for  normal 
stresses  in  the  1  direction). 

Y2  axial  flow  stress  along  the  transverse  material  direction  of  a 

transversely  isotropic  material  (for  normal  stresses  in  the  2  and  3 
directions  ). 

Y4  shear  flow  stress  in  the  isotropic  (i.e.  transverse-transverse )  plane 

of  a  transversely  isotropic  material  (i.e.  for  shear  stresses  in  the  2-3 
plane);  contracted  form  of  Y23. 

Y3  shear  flow  stress  in  a  plane  normal  to  the  isotropic  plane  of  a 

transversely  isotropic  material,  known  as  the  longitudinal-transverse 
shear  strength  (i.e.  for  shear  stresses  in  the  1-2  and  i-3  planes); 
contracted  form  of  Y13. 

Ej  Youngs  modulus  in  direction  i. 

Shear  modulus  in  i-j  plane. 

v’1j  Poisson’s  ratio  in  i-j  plane. 
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